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ABSTRACT 



Experiments aimed at detecting highly-redshifted 21 centimeter emission from the 
Epoch of Reionization (EoR) are plagued by the contamination of foreground emission. 
A potentially important source of contaminating foregrounds may be Faraday-rotated, 
polarized emission, which leaks into the estimate of the intrinsically unpolarized EoR 
signal. While these foregrounds' intrinsic polarization may not be problematic, the 
spectral structure introduced by the Faraday rotation could be. To better understand 
and characterize these effects, we present a simulation of the polarized sky between 
120 and 180 MHz. We compute a single visibility, and estimate the three-dimensional 
power spectrum from that visibility using the delay spectrum approach presented in 
Parsons et al.l (l2012bl ). Using the Donald C. Backer Precision Array to Probe the Epoch 
of Reionization (PAPER) as an example instrument, we show the expected leakage into 
the unpolarized power spectrum to be several orders of magnitude above the expected 
21cm EoR signal. 



Subject headings: Cosmology: observations, Instrumentation: interferometers, polar- 
ization 



1. Introduction 



A significant amount of thought has gone into the problem of foreground removal allowing for 
the detection of the power spectrum of neutral hydrogen during the Epoch of Reionization (EoR) 



if 



(|Bowman et al.ll2009l : iMorales et al.ll200fil : iLiu et al.ll2009l : iLiu fc Tegmarkll201ll : iDillon et al.ll2012l . 
e.g.). Essential to nearly all these techniques is the spectral smoothness of these foregrounds, which 
allows for a separation in k space of the foreground emission and the signal from the EoR. Various 
mechanisms will leak polarized emission into the best estimate of the EoR signal, which then may 
introduce spectral structure to an EoR experiment's measurement. 
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By contrast to unpolarized foregrounds, co mparatively little work has addressed the problem of 
detecting and removing polarized foregrounds. I Pen et al.l (|2009l ) provided some of the first relevant 
upper limits at small angular scales using the Giant Metre- Wave Telescope (GMRT). They found, 
for spherical harmonic multipoles 200 < I < 5000, an a ngula r power spectrum upper limit of around 
Ce < 100 mK 2 . More recent work by iBernardi et al.l ( 20101 ) detected polarized power at the same 
level at £ < 1000 using th e Westerbork Synthes is Radio Telescope (WSRT), with no significant 
detection above I of 1000. 
polarized point sources. 



Pen et al 



Bernardi et al.l (120101 ) did not detect emission directly attributable to 



(2009) present their upper limit in terms of the three-dimensional 



power spectrum, but IBernardi et al.l (|2010l ) only present a C% spectrum which has been integrated 
along the frequency direction. Figure Q] gives a summary of the low-frequency measurements of 
polarized power spectra. It is clear that the angular power spectrum of the unpolarized sky mus t 
be scaled by a mean polarization fraction of about 0.3% to agree with the IBernardi et al.l (J2010) 
measurements. This requires a significant degree of depolarization of synchrotron emission from 
ordered magnetic fields. 



Jelic et al.l (J201Cj) attempted to further constrain the problem by fully simulating a full-Stokes 
realization of galactic synchrotron emission over a 10° x 10° field of view. They present a realistic 
spectrum of the mean temperature of polarized emission, but do not extend their analysis into 
the power spectrum. Their analysis also predicts a polarization frac tion from diffuse emission 
much higher than the limited measurements ava ilable. iGeil et al J (1201 If ) also investigate the issue, 
proposing the use of the RMCLEAN algorithm (jHeald et al.l 120091 ) to mitigate the effects of these 
polarized foregrounds. While these papers provide detailed descriptions of the polarized sky and 
removal strategies in the image plane, they provide little discussion of the line-of-sight direction. 

This paper aims to steer the discussion of polarized foregrounds towards the terra incognita of 
the third, frequency dimension. We begin in section [2] reviewing the basics of polarize d interferome 



try, and the delay spectrum approach to estimating the power spectrum, presented in lParsons et al 



(|2012bl ). We discuss the relevance of polarized foregrounds to measuring the 21cm power spectrum. 
In section [3l we discuss the design and implementation of a simulation of the polarized sky and its 
results. Finally, in section [4] we briefly discuss prospects of polarized source removal and leakage 
mitigation. 

Throughout the pape r, we will use the Do nald C. Backer Precision Array to Probe the Epoch 
of Reionization (PAPER) (jParsons et alj|2010l ) as a model instrument. The results are not spe- 
cific to that instrument, nor are they specific to the Delay Spectrum analysis used in this pa- 
per (jParsons et al.l l2012bl ). Any 21cm EoR p ower spectrum det ection experiment with linear 
feeds, including; t he Murchison Widefield Array (ITingay et al.ll2012l ) or the Low Frequency ARray 
( Rottgeringlhood ) could fall subject to the leakages described here without a perfect calibration. 
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2. Preliminaries 

2.1. Definition of Stokes Visibilities 

Two prominent ways in which polarized sky emission can leak into an interferometeric esti- 
mate of Stokes I are leakages due to non-orthogonal and rotated feeds and beam ellipticity — an 
asymmetry in the two linear polarizations of a primary beam which causes unpolarized signals to 
appear polarized, and vice versa. Th e first is a well-understood question, discussed in length in 



the series of papers bv lHamaker et aL| (j 19961 ). This type of leakage can be corrected by the proper 
linear combination of visibilities. Hence, we will focus on the latter issue. To begin, we will examine 
the contents of an interferometric spectrum, and relate them to the intrinsic Stokes parameters. 

As a reminder, we present the measurement equation for an interferometer in the flat-sky limit: 
V ab (u,v,u) = [ A ab (l,m)I ab (l,m)e- 2 ^ uv+lm ^ dldm. (1) 



Here the polarization indices a, b indicate the polarization state of the measurement. A ab is the 
primary beam, I and m are direction cosines of the celestial sphere (with their Fourier components 
u and v), and I ab is the sky emission projected along the polarization state of the measurement. 
Henceforth, we will be writing all visitibility equations in the flat-sky limit. We can do this without 
loss of generality since the polarization properties of a visibility are unaffected by this assumption. 

It is worth noting that each linear feed measures a one-dimensional projection of the incident 
electric field. This causes correlations with different feed orientations to contain information about 
different polarization states of the incident radiation. A convenient short-hand notation for the 
polarization content of each visibility is 



~^xx Vxy \ 

VyX VyyJ 




(2) 



The Jones matrices Jjj (see lHamaker et al.lll996l ) relate the sky emission (/, Q, U, and V) to an 



interferometeric measurement. Each antenna's Jones matrix is dependent on nearly all instrumental 
parameters, but for the purposes of this paper, we will investigate the effects of direction-dependent 
gains, or the primary beam. 

In the flat-sky approximation, an interferometer natively measures the two-dimensional spatial 
Fourier transform of the sky. Ideally, this would allow an observer to estimate the three-dimensional 
power spectrum by simply transforming the frequency, line of sight direction, and cross-multiplying 
measurements without imaging. Relaxing the imaging requirement provides an incentive to make 
estimates of the Stokes parameters — / for instance — in the visibility domain, where the Stokes 
parameters are not defined. A naive addition between the linearly polarized, xx and yy visibilities 
should estimate the total power of sky emission, Stokes /. Similar operations can be performed for 
all polarization states. A sensible method of estimating the four Stokes parameters in the visibility 
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domain is to add visibilities as images are typically added. Hence, we define Stokes visibilities (for 
linear feeds) as: 
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Explicitly writing the expression of V/ by substituting equation [T] into equation [3l we find that 

Vl =VxX + Vyy 

= J (A xx + A yy ) Ie~ 2 ^ lu+vm) dldm 

+ J (A xx - Ayy) Q e -^(lu+vm) dldm (4) 

Writing the visibility this way highlights one source of polarized leakage: that due to beam ellip- 
ticity. If the xx and yy beams are not equal, the last term of the equation [5] would be non-zero. 
Polarized emission then enters the estimate of /, weighted by the differenced beam. This equa- 
tion also points out the difference between Stokes parameters and the Stokes visibilities defined in 
equation [3l They estimate their respective Stokes parameters, but still include terms from other 
polarization states. 

Thus, if an instrument's linearly polarized beam is not symmetric under a 90° rotation, all the 
intrinsic linear polarization will not cancel, corrupting the visibility estimate of Stokes parameters. 
Figure [2] shows two slices through the model beam for PAPER, through constant azimuth and 
through zenith, and a cut through the same beam, rotated by 90°. We note that these slices are 
not identical, and show the differenced beam, which provides the mechanism for leakage. 

In an image-based analysis, these effects may be corrected by a simple reweighting of the image, 
but these reconstructions are subject to the accuracy of the primary beam model. The results of 
this paper indicate the leakages due to the primary beam shown in figure [21 but could equally be 
applied to an image-based analysis whose primary beam is only known to the level of the differences 
in that figure. 

Once more, we note that the Stokes visiblities are only the best guess at the true Stokes 
parameters. Without an exact characterization of the primary beam or dense u-v sampling, one 
does not have the ability to fully correct for the leakages mentioned. Hence, it is imperative to 
inspect the Stokes I signal's corruption. 
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2.2. The Delay Spectrum Approach 



The delay-spectrum approach to measuring the 21cm power spectrum ([Parsons et al.ll2012bl ) 
requires no imaging, again providing an incentive for adding raw visibilities. The delay-spectrum 
approach embraces the natural units of an interferometer by sampling the power spectrum in base- 
line and frequency dimensions, native to an interferometer. With this approach, each baseline can 
be individually transformed into an estimator of the power spectrum of the incident temperature. 
This method also prevents the small-s cale structure introduced by Fourier transforming over gain 



calibration errors ([Morales et al.ll2012l ). 



Perhaps the greatest advantage of the delay spectrum approach is its ability to isolate smooth- 
spectrum foregrounds on a per-baseline basis. This relaxes the burden of isolating sources in the 
image plane, and all ows for a more sparsely sampled array (e.g. a redundantly sampled one, as in 



Parsons et al.ll2012al ). The 21cm EoR's extent to super-horizon delays depends precicely upon its 
spectral non-smoothness. It is imperative, then, to identify and characterize spectrally non-smooth 
foregrounds that may corrupt the EoR measurement. 

The delay spectrum approach is particularly useful for simulating the effects of a systematic 
error on the power spectrum. Because each visibility estimates the power spectrum, a simulation 
needs only create one visibility, which encapsulates much of the power spectral information needed 
for analysis. The computational cost is additionally reduced by the ability to calculate a power 
spectrum by Fourier transforming only one visibility, rather than gridding several visibilities into 
a (u, v, v) cube and Fourier transforming along each u, v pixel. 

It is important to note that the results of this paper are not limited to delay-spectrum-based 
analyses. All effects mentioned in this paper are native to the sky and its expected spectral 
structure. Any analysis technique will incur all the same issues; we use the delay spectrum as a 
convenient probe for these issues. 

We briefly review the definition of the delay transform as a Fourier transform over the frequency 
axis of a visibility: 

V(r)= / V{u,v,v)e~ 2 ™ TV dv, (5) 

where the tilde denotes a delay-transformed visibility. We treat this as an estimator of the 
Fourier transform of the temperature field, and squaring it, noting that 

V 2 (T) = |V(T)| 2 ocP(fc). (6) 



To better mimic actual measurements, we divide the band into ten, 6-MHz sub-bands, and 
perform the delay transform on each individual band. The sub-bandwidth is chosen to sample 
the maximum cosmological distance over which the 21cm signal is expected to be cotemporal 
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(jWvithe fc Loebll2004l ; iFurlanetto et al.ll2006al . e.g.). Once these transforms are computed, we mul- 
tiply each spectrum by the appropriate factors to obtain the power spectrum in units of temperature 
squared. We write the expression converting the squared, delay-transformed visibility V(r) into a 
"unitless," cosmological power spectrum A 2 (k) as 



V 2 (r) 



QB 2vr 2 



2k B \ 

X 2 J X 2 Y k 3 



A 2 (k), 



(7) 



where ks is the Boltzmann constant, A is the observing wavelength, £1 is the solid angle of the 
primary beam, B is the observed bandwidth, and X and Y are cosmological scalars which convert 
observed angles and frequencies in to /iMpc -1 . 



Parsons et al.l (|2012bl ) offers a much more detailed discussion of this approach, beyond the 



scope of this paper. 



2.3. Sparse u-v Sampling and Wide-Field Polarimetry 

Another advantage of the delay spectrum approach is that it relaxes the requirement of gridding 
in the u-v plane. Each baseline is assigned a position in the u-v plane ab initio, and visibilities from 
similar baselines may be coherently added without imaging. This allows for sparse sampling in the 
u-v plane without damaging effects from side lobes or missing data, problems other methods may 
experience. Since the delay spectrum rotates a power-spectrum estimate int o the native coordinate 



system of an interferometer, there are no inherently missing frequency-data. iParsons et al.l (|2012al ) 
presents the sensitivity benefits of a sparse, redundant array configuration, but other techniques 
aim to uniformly sample the (u, v, v) cube, in order to mitigate the systematic effects of computing 
a Fourier Transform across unevenly sampled data. 

An obvious disadvantage of having sparsely-sampled data is poor imaging. Not only does sparse 
sampling provide a highly irregular synthesized beam, but it also limits the available information 
for a full reconstruction of the image. Without adjacent u-v samples, a full, accurate deconvolution 
by a wide beam simply has insufficient information. As we will see, the inability to correct for 
beam effects will provide a significant source of systematic error via polarized leakage. 

By choosing to wield the full power of the delay spectrum approach and redundant sampling, 
an observer is forced to add visibilities with no beam weighting. The beam information supplied 
by adjacent u-v samples simply does not exist, and without transforming into the image plane is 
unrecoverable. Hence, the imperative to investigate the implications of a lack of beam-weighting 
arises. 

Together, redundant sampling and the delay spectrum approach give a 21cm EoR experiment 
incentive to add raw visibilities, subjecting it to potential leakage. An elliptical primary beam 
gives a mechanism whereby polarized emission can corrupt an estimate of the total power. To 
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what degree does polarized emission corrupt an estimate of the 21cm EoR signal. We will start by 
characterizing any spectral non-smoothness that will possibly arise from polarized leakages. 



2.4. Faraday Rotation 



Faraday Rotation affects the polarization properties of an electroma gnetic wave traveling 
through a plasma containing a magnetic field (jRvbicki fc Lightmanl Il979l ). The circular polar- 
ization oriented in a right-handed fashion to the direction of the incident field will be slowed by 
the plasma. This causes a rotation of the E-field's polarization angle by 



= A 2 . 6 0N0 f B n (s)n e (s)ds = \ 2 RM. (8) 
(m e c z y J 

Equation [8] defines the rotation measure (RM) by which we characterize this phase wrapping. 
Since the Stokes parameters characterize the square of the electric field, the phase of the polarization 
vector is shifted by twice the angle defined in equation [8l After a signal passes through a Faraday 
screen, we measure a rotated polarization angle, 

(Q + iU) meas = (Q + iU) mt e-^ 2 , (9) 

where the subscript "meas" denotes the measured Q and U measurements, and "int" denotes the 
intrinsic Q and U signal, as would be measured behind the Faraday screen. 



3. Simulation 

3.1. Single-Source Power Spectrum 

We begin our investigation of the effects of polarized foregrounds on the 21cm EoR signature by 
examining the power spectrum of a single source at zenith, whose signal has the structure of a single 
RM. In doing so, we can develop an intuition for the rotation measures that affect cosmologically 
interesting k modes of the power spectrum. By looking at what is effectively the impulse response 
of a Faraday screen on the power spectrum, it will be easier later to interpret a more complicated 
model. 

Figure [3] shows the spectra of a few linearly polarized sources behind Faraday screens, S(u) = 
exp{— 2iRM\ 2 }. Each spectrum contains one, one-Jansky source, located at zenith (delay of zero). 
Notice that at the highest rotation measure shown, the spectrum is not critically sampled at the 
lower frequencies. This is due to the uneven sampling of A 2 across the band: as AA 2 ~ dX 2 /dv Av ex 
Avjv^ increases, the sensitivity to large rotation measures decreases. 

Figure 0] shows the delay-spectrum over PAPER'S full 100-MHz band of the same Faraday 
screens as Figure HI Each of these spectra samples several cosmologically disconnected regions, 
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but we note that these measurements are not cosmologically significant, only revelatory of Faraday 
leakage. The results of these transforms over the full simulated band are shown in figure [5j The most 
important feature of this plot is this: there is a k mode associated with each rotation measure at 
each redshift. We can construct an analytic estimate of this k mode by setting the argument of the 
exponents of a delay mode and a rotation measure mode to sum to zero. First, we approximate the 
cosmological /c-mode sampled by a delay mode as r ~ kudru/du. Next, we recall the cosmological 
scaling from frequency into /iMpc -1 , 

dr\\ _ dr^dz_ _ c(l + z) 
dv dz dv H(z)v 

Finally, we set h\ ■ dru/du ■ v + 2RMX 2 = 0. Substituting Equation \W\ for the derivative, we derive 
an expression for the k-mode most affected by a rotation measure RM: 

k p^.MA 2 , (11) 
1 cl+z v ' 

where H{z) is the Hubble parameter. Figure [6] shows a plot of the most culpable rotation measure 
versus frequency and redshift. 



3.2. Full-Sky Simulation 

To better grasp the effects of Faraday leakage into the 21cm signal, we generate several random 
realizations of the sky, each consisting of many polarized point sources. Each source passes through 
a Faraday screen with a randomly assigned rotation measure. Next, we simulate that source for 
a single baseline. Finally, we calculate the power spectrum measured by that visibility. Only one 
visibility needs to be simulated, because the delay-spectrum approach makes use of the fact that 
each baseline measures the 21cm EoR with a range of /c-modes determined by the baseline length, 
orientation, and bandwidth. 

As mentioned before, we use PAPER as a model for choosing instrumental parameters. Our 
simulated band goes from 120 to 180 MHz, with 60 kHz channels. We use a simulated primary 
beam model of the PAPER dipole, which has a full width at half max of around 45° at 150 MHz 



(jPober et alJl2012h . 



Rather than creating an exact simulation of the physical sky, we create a simulation whose 
statistical properties are physically motivated. This choice reflects a desire for simple, easily tunable 
parameters for the simulation. In that same spirit, we model all sources simply as point sources with 
a Poisson distribution. The simulation's primary concern with the spectral information of polarized 
foregrounds allows us to justify neglecting the angular terms. This is equivalent to assuming for all 
the relevant k modes, fen S> k±. Emphasis on the fe|| or spectral modes also motivates our decision 
to model the sky a s numerous point sources. For a more detailed discussion of these effects, we 



direct the reader to I Jelic et al.l (|2010i ) . 
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Source positions are distributed uniformly over the sphere. A single source's altitude 9 is 
drawn from a distribution in which cos 9 is uniform on [0, 1]. A source's azimuthal angle 4> is drawn 
independently from cos 9 from a distribution uniform on [0, 2ir]. This choice of source position 
distributions conserves the density per area of sources across the sky, and is equivalent to drawing 
both direction cosines, (l,m) = (sin 9 cos <f), sin 9 sin (f>) , from a uniform distribution on [—1, 1]. 

In order to achieve realistic source fluxes and source counts, we base the distributions from 
which we dr aw various parame ters on previous radio surveys. For t he source fluxes, we aim to agree 



with VLSS (jCohen et al.ll2007h . NVSS surveys (|Condon et al.lll998h . and the 6C survey (Hales et al 



1988). For the polarization information, we aim to agree with polarized measurements taken by 
the NVSS survey, particularly, their nearly full sky rotation measure map (jTavlor et al.ll2009l ). as 
well as the polarized survey of the VLSS. 

We present t wo scenarios for the source counts. First, we draw from the 6C distribution 
(IHales et al.lll988Tl. taken in the PAPER band at 151MHz. Second, we extrapolate VLSS source 
counts (ICohen et al.l 20071) from 74 MHz to 150 MHz, using a spectral index of -0.79, following the 
work of I Cohen et al.l (|2004l ). The latter source counts provide more sources at higher flux, which 
infect the power spectrum from the / visibilities, as we will show in the following sections. Both 
number counts are c onsistent with the recent measurements by the Murchison Widefield Array 
([Williams et al.ll2012l ). an instrument similar in many regards to PAPER. 



The source counts found bv lHales et al.l (|1988l ) may be characterized by two power laws, turning 
over at some knee flux S 



dN_ _ J4000 S-°- 76 S- L75 Jy^sr- 1 S min < S < S a 
dS ~ [4000 S~ 2m Jy^sr- 1 S Q < S 



(12) 



Following the 6C survey, we choose the turning point, S Q to be 0.88 Jy. The number of sources 
simulated (14,855) is chosen by the size of the PAPER beam at 151 MHz (0.76 sr) and a flux range 
over which to integrate. We choose to include those sources in between 100 mJy and 10 Jy. This 
choice provides a reasonable dynamic range of sources. Below the lower limit, the 6C sources are 
unreliable, and we assume sources above 10 Jy may be easily identified and removed. 

If we were to blindly extrapolate the 6C source counts below the lower limit of the catalog, we 
would add a negligible amount of power. Integrating S 2 dN/dS down to some minimum flux esti- 
mates the contribution of the sources above that flux to the total variance of flux. By inserting the 
6C source counts, we find that we are including ~ 70% of the estimated total variance. Extending 
the minimum flux would add more power to the simulation, but would not drastically alter the 
results of this paper. 



VLSS source counts follow a single power law, given bv lCohen et al.l (|2004l ) as 

4865 S~ 2 - 3 Jy- 



dN 
~dS 



■-\sr-\ 



(13) 
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For these source counts, we choose minimum and maximum fluxes of 0.8 Jy and 100 Jy, 
respectively, rejecting sources well below the lower limit of the catalogue, and providing a reasonable 
dynamic range for the included sources. Integrating over the PAPER beam provides 6995 sources. 
These source counts are not qualitatively different from the 6C counts at low fluxes, but they do 
differ substantially in that they provide many more bright sources. 

For both source counts, we also assign a spectral index to each individual source spectrum, 
drawn from a normal distribution, mean -0.8, standard deviation 0.1. 

Each source is also assigned a random polarization angle %i uniformly s ampled from [0, ir 



Total flux is multiplied by a polarization fraction, chosen to reflect the studies of lTucci Toffolatti 



([2013). We sample the polarized fraction (II) from a log-normal distribution whose mean is 2.01% 
and whose standard deviation is 3.84%. Because the log-normal distribution is not upper-bounded, 
we reject any drawing over 30%. Following the aforementioned study, we do not impose any corre- 
lation between source flux and polarization fraction. It has been noted that, among ot her effects 



band width depolarization causes the polarize d fraction to decr ease at lower frequencies ([Law et al 
201 ll ). This along with measurements from |Pen et al.l (|2009l ) indicates that these distributions, 



taken at 1.4 GHz, may overestimate the distribution at 150 MHz. We neglect these effects, taking 
the 1.4 GHz distribution at face value, since the mean polarization fraction can be thought of as a 
scale factor to the overall power spectrum. 

The distribution of Rotation Measures for these sources influences the power spectrum perhaps 
more than any other parameter. Unfort unately, assigning a physically motivated distribution to 
these is not such a straight-forward task. I Taylor et al.l (120091) generate a map o f Rotation Measures 



for all lines of sight above declination —40°. However, as lPshirkov et al.l (|201ll ) mention, each line 



of sight contains a significant amount of un certainty. They compare Taylor with the more recent 



study by lKronberg &: Newton- McGed (|201lh and find large disagreements. Furthermore. lLaw et al 



(|201ll ) point out the depolarizing effects of coarse spatial resolution and a wide bandwidth. They 
conclude that the results of rotation measure synthesis (RMS) is highly instrument-dependent. Low 
spatial resolution will average over many pointings of a potentially turbulent magnetic field, effec- 
tively averaging the measured rotation measure down. They also mention two effects of bandwidth 
depolarization. First, they note that at low frequencies, any reasonable spectral resolution will sum 
many wrappings of a rotation measure component, reducing the measured power, and more highly 
resolved spectrometers tend to measure higher rotation measures. 

Taking these effects into consideration, we use a distribution based on the NVSS Rotation 
Measure maps, smoothed to 8° scales. This averages over much of the uncertainty, and (hopefully) 
generates a distribution more similar to what exists. Because PAPER has at the center of its band 
a ~ 20' synthesized beam, we will neglect the effect of higher rotation measures, like those measured 
in Kronberg and Newton-Mcgee. The width of PAPER'S frequency channels causes the instrument 
to be insensitive to such high rotation measures, and their effect would be noise-like. 



Histograms of the distributions of rotation measure, polarized fraction, and source counts can 
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be found in figure Over-plotted on all is the distribution from which they are drawn. 

We model all sources as point sources, neglecting the effects of any diffuse emission. This choice 
reflects the desire for a simple, easily tunable parameters in the simulation. While diffuse emission 
certainly is present, its spectral structure is qualitatively the same as that of a point source, so the 
results in the frequency direction will not change. Were we to add a diffuse model to the simulation, 
it would widen a source's response in delay-space, preventing it from being localized. This will not 
affect the line-of-sight power spectrum, at high k, since the spectrally-smooth, diffuse emission still 
falls within the horizon. 

Table [T] summarizes the three treatments of the simulation we will be using. Simulation A, 
with 6C source counts and the Taylor RM distribution is likely the most accurate. Simulation B 
steepens the source counts, providing fewer, brighter sources, and simulation C doubles the width 
of the RM distribution. 



To check if the results of this simulation are consistent with the measurements in Bernardi et al 



(|2010|), we compare the two-dimensional Cg power spectrum with that presented in their paper. 



Figure [8] shows t he 2D power sp e ctrum of one realization of the simulation, with the constraints 



from figure 20 of iBernardi et al.1 (|2010l ) plotted over it. We see qualitatively that our simulation 



well obeys the upper-limit imposed, and proceed with the results. 

We calculate the visibilities for a 32m, east-west baseline. T his choice reflects the m ost common 



spacing of the maximum-redundancy configuration presented in iParsons et al.l (|2012al ) . The choice 
of baseline orientation is arbitrary, and since we are modeling only point sources, the choice of 
baseline length will only set the horizon limit of the power spectrum. Since the delay affected by 
a rotation measure is independent of a choice of baseline (equation II lj) . choosing a relatively short 
baseline will isolate the foregrounds at lower r, and highlight the Faraday leakage. 

The full measurement equation used in this simulation is 



'ISOMHzV 



V(u,v,u) = ^2A(l J ,m„u)U j S} 50 l 

i=i \ v J 



x exp { -i[2irv(ulj + vrrij) + 2RMjX 2 + 2xj] } , (14) 

where each source j is assigned a flux (Sj), polarization fraction (II?), spectral index (ay), a position 
(lj,rrij). rotation measure (RM), polarization angle (xj)> an d is weighted by the model primary 
beam (A). A sample Q visibility is shown in Figure [9l 



Name 


Source Counts 


A'src 


RM distribution 


A 


6C 


14,855 


NVSS 


B 


NVSS 


6,995 


NVSS 


C 


6C 


14,855 


2 x NVSS 



Table 1: Simulation Treatments 
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We choose not to include the parallactic rotation of Q into U, implying that the Q we label in 
this paper are fixed to topocentric, azimuth and altitude coordinates. This choice clarifies equations 
and allows for an ease of understanding which would be obfuscated by writing both Q and U. 



3.3. Simulated Power Spectra 

Figure [10] show the power spectra for source counts from the extrapolated VLSS and 6C 
surveys, as well as the spectrum of 6C source counts with a widened rotation measure distribution. 
We interpret the power spectrum of the I visibility as the amount of polarized leakage corrupting 
the EoR signal (henceforth called Q —■ I leakage), and the Q visibility's power spectrum is our best 
representation of the polarized signal. These plots show the median power in each k bin for 1000 
realizations of the simulation, with error bars show the one-sigma extent of the bandpowers for these 
realizations. These power spectra confir m the prediction of sec tion 13.11 that A 2 phase wrapping 



extends the foreground cutoff presented in lParsons et al.l (|2012bl ) to higher r bins, corrupting some 



of the most sensitive regions of k space for 21cm EoR analysis. They also demonstrate the prediction 
in that section that high-redshift bins will be most affected. 

The severity of the leakage can be inferred from the power in the most EoR-sensitive k bins 
(0.2 hMpc- 1 <k<0.3 /iMpc" 1 ). Figure [TT1 shows A 2 (fc) in these bins as a function of redshift. 
The leaked power ranges in the thousands of mK 2 , increasing from high frequency / low redshift to 
low frequency / high r edshift. These es timates are about two orders of magnitude above level of the 



expected 21cm signal (jLidz et al.ll2008l ). If we take this simulation as an accurate prediction of the 
low- frequency sky's polarized emission, these results imply that naively adding V xx and V yy , formed 
with an approximately 10% asymmetric primary beam, incorporates enough bias from polarized 
leakage to completely obscure the 21cm signal. The levels of leakage in our simulations demands a 
strategy to model and remove polarized sources. 

We note that simply adding V xx and V yy will also remove a negligible component of the EoR 
signal via the same mechanism. In a sense, Q — > I leakage can be thought of as a rotation of 
power between the two Stokes parameters. Hence, for precision measurements of the EoR signal, 
this simple estimate may not be ideal. However, the effect of I — > Q leakage is small (compare the 
levels of Vq and high-A; modes of V/) and should not provide a significant hinderance to detection. 



4. Mitigation 

The previous section predicts an excess polarized signal due only to point sources at around 
10 4 mK 2 at k ~ 0.15 for most treatments of the simulation. While the exact levels of these 
predictions may be subject to some error, the need certainly arises for some removal scheme, which 
must roughly suppress power from polarized foregrounds at around the -40 dB level in the power 
spectrum. 
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To investigate the effects of modelling and removing polarized sources, we rerun the simulation, 
excluding the brightest polarized sources. Figure [12] shows the median value of 500 simulations of 
one /c-bin of the power spectrum for removing the brightest 1000, 2000, 3000, and 5000 sources. 
These limits in number of sources correspond to flux-limits of 40 mJy, 25 mJy, 18 mJy, and 11 mJy, 
respectively. We remove these sources from treatment A of the simulation, which includes around 
fifteen thousand sources. Despite having removed nearly one-third of the sources, the leaked power 
still exceeds 10 mK 2 . 

To remove enough flux to consistently fall below the expected EoR signal, we needed to remove 
eight to ten thousand sources (~6 mJy), roughly the number of sources in the simulation. We 
recomputed these fluxes with a lower minimum flux (60 mJy), expecting a similar result, but found 
that we increased the power in this fc-bin by one or two mK 2 . Ten to twelve thousand sources 
required removal for the total power to fall below 10 mK 2 . We exclude further investigation of 
this analysis for three reasons. First, current measurements do not constrain dN/dS to the levels 
necessary to accurately model such low-flux sources. Second, including lower-flux sources does not 
significantly affect the result that the expected polarized power spectrum will be of the order of 
10 4 -10 6 mK 2 . Third, the variance from one simulation to the next in the power was large enough 
that the two treatments of the simulation — even with ten thousand sources removed — could not 
be considered significantly different. 

These onerous levels of source-removal suggest that a different mitigation scheme be considered. 
Future instruments may take polarization into consideration in design. Leakage can be mitigated 
with more circular beams, and c ircular feeds avoid Q — >■ J leakage entirely. Even with existing 



data, rotation measure synthesis (jBrentjens h de Bruvnll2005l ) could potentially provide the power 



to separate sources with distinct RM structure to be separated from EoR signal. 



5. Discussion 

Were a power spectrum computed from only one linearly polarized visibility (xx, for instance), 
all polarized power would corrupt the measurement. We have chosen to suppress the polarized 
leakage by adding linearly-polarized visibilities. The leakage is dependent on the difference of the 
two beams, and by having beams that are at most ten percent different suppresses the signal by 
around -30 dB. Correcting for the beam-weighting in the image domain can further suppress the 
leakage, but errors in the beam model will introduce leakage in much the same manner. Hence, 
the constraint of having to suppress polarized leakage by -40 dB causes the need for an accurate 
primary beam model to around the one percent level in the case of imaging, or symmetric at the 
one percent level if visibilities are used directly. 

These estimates of power are also dependent on the relative strengths of diffuse, polarized 
emission and polarized point sources. He have taken care to agree with current measurements, but 
we note that above £ ~ 300, the current constraints are noisy. We have interpreted them as an 



-14- 



upper limit. In the limiting case where diffuse emission is the only component to the polarized sky, 
this leakage could be suppressed by measuring with a longer baseline, which in turn measures lower 
i and k±. We have chosen to use a sixteen- wavelength baseline, which corresponds to t ~ 200. 
This choice of baselin e length is relatiyely s hort for interferometers at these wavelengths, but falls 
at the high end of the iBernardi et al.1 (|2010l ) detection. 



Including diffuse emission in the simulation would slightly alter the amplitude of the polarized 
power spectrum, although it would not eliminate the frequency structure. The greatest effect of 
including diffuse e mission would be th at it selects a range of rotation measures to include in the 
measurement. ThelTaylor et al.l (120091) map shows corre lation of rotation measures on large scales 



(jPshirkov et al.l l201ll ; iKronberg fc Newton-McGed l201ll ) , which would allow only a certain range 



of rotation measures to enter a measurement with each pointing. Choosing a particularly benign 
region of the polarized sky could potentially mitigate some of this leakage. The footprint of diffuse 
and point-like emission in the frequency direction is identical - diffuse emission would only affect 
the k± dimension of the power spectrum. 

We conclude our discussion of the simulation results by noting the large variance in the sim- 
ulated power. The results shown are the median band-powers in A 2 (fc) for 500 realizations of the 
simulation. Taking so many realizations into account essentially maps out the posterior distribu- 
tion of the A 2 (fc) bandpowers. The one-sigma width covers nearly an order of magnitude, which 
indicates the level of Faraday leakage is highly sensitive to the exact parameters drawn in any 
realization. The actual level of leakage measured will thus be highly dependent on a choice of field, 
and on cosmic variance. 



6. Conclusion 

We predict the three-dimensional power spectrum of polarized emission around 150 MHz to 
be in the range of 10 4 -10 6 mK 2 at k\\ ~ 0.15 /iMpc" 1 . These predictions are based on simulations 
movitvated by current observations of the polarized sky at 150 MHz and 1.4 GHz. An elliptical 
beam provides one mechanism for this power to leak into a measurement of the unpolarized signal. 
Using a fiducial model of the PAPER beam, we estimate this leakage to be in the thousands of 
mK 2 , several orders of magnitude above the expected 21cm EoR signal. Modelling and removing 
polarized sources may eliminate much of this leakage, but these simulations suggest the amount of 
removal required far exceeds reasonable capabilities of current instruments. 

Work is currently underway to measure frequency structure of polarized power, to investigate 
the amount of its leakage into the I power spectrum, and to better characterize the polarized radio 
sky using existing PAPER data. 
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Fig. 1. — Recent es timates of low-frequen cy polarized power spectra. The thick black points with 



error bars show the Bernardi et al 



shows the upper limit of 



Pen et al 



1 2010) measurements of a field around 3C1 96, and the solid lin e 



((20091). The Haslam map at 408 MHz (jHaslam et alJ 11982 ). 
scaled by a polarization fraction of 0.3%, is shown by thin points, and a power-law extrapolation 
is shown with a dotted line above I = 200. This fraction was chosen to agree with the low-£ points 
in the Bernardi measurement. At high-^, the upper limits do not constrain the level of polarized 
emission. The grey vertical line shows the £ mode sampled by the simulation in this paper, as a 
point of reference. 
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Fig. 2. — Power vs. altitude through zenith for east-west slices of the xx and yy PAPER beams 
(dashed, dotted respectively). The solid line shows their difference. 
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Fig. 3. — Real component of the visibility V = exp{— 2RMX 2 } for rotation measures 3, 10, 30 and 
100 m~ 2 , from top to bottom. This plot is intended to demonstrate both the frequency-dependent 
phase wrapping and the frequency-dependent critical sampling of faraday rotation. Notice how the 
100 m~ 2 component (located at the bottom) is critically sampled at high frequencies, but loses 
coherence at low frequencies. 
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Fig. 4. — Delay spectra of the Faraday-screened visibilities shown in figure El This plot demon- 
strates the effect of a Faraday screen on the delay spectrum of a source. There are two effects: 
first, a faraday screen widens the response of a spectrum in delay space, and second, the faraday 
screen scatters power to high delay. This causes a potential problem as many foreground removal 
techniques require smooth-spectrum foregrounds. Note that higher rotation measures mimic higher 
delays. 
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Fig. 5. — Power spectra of the visibilities plotted in Figure El computed for the median redshift bin 
of the PAPER band (z ~ 8). This plot demonstrates the coupling between rotation measure and 
fc-modes, and confirms that Faraday-rotated spectra scatter power to higher k modes than would be 
contained within the horizon. Grey, vertical lines show the maximally-infected /c-mode, predicted 
by equation [TTJ The rise in power at high k is due to the ~ 10 -9 sidelobe of the Blackman-Harris 
filter. The units of the y-axis are tuned so the total power integrates to one. 
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Fig. 6. — The rotation measure most affecting k ~ 0.15 /iMpc , calculated from equation 
Shaded region indicates the range of rotation measures affecting 0.1 < k < 0.2 /iMpc~ . 
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Fig. 7. — Distributions of simulated parameters. (Top Left) S 3 / 2 N(> S) source counts of unpo- 
lariz ed flux. In black, th e 6C source counts, from lHales et al. ( 19881 ). and in grey, the VLSS counts, 



from lCohen et al.1 (120071 ). ( Top Right) Norma li zed h istogram of the log-normal distribution of po- 
larized fraction, taken from lTucci &: Toffolattil (120121 ) . (B ottom) Norm a lized histograms of the two 
distributions of rotation measures. Black is taken from iTaylor et al.1 (120091 ) . grey is that with a 
doubled standard deviation. In each plot, the dashed line is the distribution from which sources 
are drawn, and the bins or points are the values of one realization of the simulation. 



-23- 



10" 



° 10 1 



10 L 




I vs . Cp f or one realization of the simulation (solid line), alongside the results from 



Bernardi et al.l (|2010l ). (crosses). This roughly demonstrates the agreement of the simulation pre- 
sented in section 13.21 with current observations. The simulation is consistently lower than the 
observed values due to the lack of extended structure. 
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Fig. 9. — Simulated Q visibility with the parameters shown in figure [71 
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Fig. 10. — Power spectral measurements for the three treatments of the simulation shown in Table 
[TJ The left column shows the power spectra of the I visibilities, and the right shows the spectra of 
Q visiblities. The top row shows simulation treatment A, the second row shows treatment B, and 
the third row shows treatment C. Line styles depict different redshift bins: 7.2 (solid), 8.0 (dashed), 
8.9 (dot-dashed) , and 10 .2 (dotted). Grey line give a toy model of the expected 21cm emission from 



Furlanetto et al.l (|2006bl ). Both visibilities include contributions from both the (intrinsic) Stokes I 
and Q. The simulated levels of Q emission indicate that polarized leakage into I needs to be less 
than -40 to -50 dB in mK 2 . 
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Fig. 11. — Median Q — > I leakage for the k ~ 0.16/iMpc _1 bin of the power spectrum, A 2 (k) vs. z 
for all treatments of the simulation shown in figure [TUl and defined in table HJ A (solid), B (dashed), 
C (dot-dashed). We note that with every treatment, the amount of polarized flux leaking into the 
power spectrum grows with z, confirming that leakage is worse at lower frequencies. 
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Fig. 12. — A 2 (k) vs z at a single k 0.16 foMpc -1 . Each line style shows the volume-averaged 
power after the brightest N sources were removed from the simulation. The number of removed 
sources, N is depicted by line style. Dotted: 1000 sources removed (40 mJy), dot-dashed; 2000 
removed (25 mJy), dashed: 3000 removed (18 mJy), and solid: 5000 removed (11 mJy). 
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